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Abstract 

A new symplectic variational approach is developed for modeling dissipation in kinetic equations. This approach 
yields a double bracket structure in phase space which generates kinetic equations representing coadjoint mo- 
tion under canonical transformations. The Vlasov example admits measure-valued single-particle solutions. Such 
solutions are reversible. The total entropy is a Casimir, and thus it is preserved. 
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Resume 

Une nouvelle approche est proposee pour modeliser la dissipation dans les equations cinetiques. Cette approche 
produit une structure a double crochet dans I'espace des phases qui conduit aux equations cinetiques d'une 
dynamique coadjointe apres transformations canoniques. L'exemple de Vlasov admet alors des solutions pour 
particule unique. Ces solutions sont reversibles ; I'entropie totale est un Casimir et est done preservee. 

Pour citer cet article iD.D.Holm, V. Putkaradze and C. Tronci, C. R. Acad. Sci. Paris, Ser. I XXX (2007). 



Version frangaise abregee Une nouvelle approche est proposee pour la modelisation des phcnomcnes dis- 
sipatifs dans les equations cinetiques [4,18]. Cette construction est realisee de telle sorte que la geometrie 
de la variable dynamique est preservee : en particulier, nous considerons I'equation de Vlasov [18] comme 
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modclc nature! dc conservation dans respace des phases. De plus nous introduisons une quantitc par- 
ticulicre nommce "mobilitc" , inspires par analogic avec la loi de Darcy pour systcmes continus avec 
auto-agrcgation [9,10,11]. Dans ce cas, on introduit la mobilitc comme le factcur dc proportionnalite 
entre la force agissant sur les particules et leur vitesse. Nous nous interessons ainsi a une forme dc dissi- 
pation qui pent gcncraliser les phcnomcnes d'auto-agrcgation aux systcmes cinetiqucs dans Tespacc des 
phases. 

Une telle approche produit une structure a double crochet [1] dans I'espace des phases similaire a 
celle presentee dans la littcraturc pour la modelisation dc certains systemes astrophysiques [13]. Cette 
structure genere une dynamique coadjointe reversible (Eq. 2) via Taction des transformations canoniques. 
On trouve finalement que toutes les fonctionncUes de la distribution sont de type Casimir et que I'entropie 
est preservee (Proposition 2.2). 

L'innovation dc notre approche se voit dans le role que la mobilitc pent jouer comme operation dc filtre 
(o moyenne) sur la fonction de distribution des particules. Par consequent on definit la mobilitc comme 
une fonctionncUc de la distribution des particules. Ce fait conduit a rexistence dc la solution dc particule 
unique qui n'est pas presente dans les ancienncs approches et represente le resultat principal de cet article 
(Theorcme 4.1). 



1. Introduction 



Non-linear dissipation in physical systems can modeled by the sequential application of two Poisson 
brackets, just as in magnetization dynamics [8]. A similar double bracket operation for modeling dissipa- 
tion has been proposed for the Vlasov equation. Namely, 
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where a > is a positive constant, H is the Vlasov Hamiltonian and [• , •] is the canonical Poisson bracket. 
When a 0, this equation reduces the Vlasov equation for coUisionless plasmas. For a > 0, this is the 
double bracket dissipation approach for the Vlasov-Poisson equation introduced in Kandrup [13] and 
developed in Bloch et al. [1]. This double bracket approach for introducing dissipation into the Vlasov 
equation differs from the standard Fokker-Planck linear diffusive approach [4] , which adds dissipation on 
the right hand side as the Laplace operator in the momentum coordinate Ap/. 

An interesting feature of the double bracket approach is that the resulting symmetric bracket gives rise 
to a metric tensor and an associated Riemannian (rather than symplectic) geometry for the solutions. The 
variational approach also preserves the advective nature of the evolution of Vlasov phase space density, 
by coadjoint motion under the action of the canonical transformations on phase space densities. 

As Otto [17] explained, the geometry of disspation may be understood as emerging from a variational 
principle. Here, we apply the variational approach to derive the following generalization of the double 
bracket structure in equation (1) that recovers previous cases for particular choices of modeling quantities. 
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Eq. (2) extends the double bracket operation in (1) and reduces to it when H is identical to E and 
/i(/) = a f. The form (2) of the Vlasov equation with dissipation allows for more general mobilities than 
those in [1,13,14,16]. For example, one may choose /i[/] — K * f (in which * denotes convolution in phase 
space). As in [12] the smoothing operation in the definition of /u(/) introduces a fundamental length scale 
(the filter width) into the dissipation mechanism. Smoothing has the added advantage of endowing (2) 
with the one-particle solution as its singular solution. The generalization Eq. (2) may also be justified by 
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using thermodynamic and geometric arguments [12]. In particular, this generalization extends the classic 
Darcy's law (velocity being proportional to force) to allow the corresponding modeling at the microscopic 
statistical level. 



2. Dissipation for kinetic equations 

We aim to model dissipation in Vlasov kinetic systems through a suitable generalization of Darcy's law. 
Indeed, we believe that the basic ideas of Darcy's Law in configuration space can be transferred to a phase 
space treatment giving rise to the kinetic description of self-organizing coUisionless multiparticle systems. 
In what follows, we will construct kinetic equations for geometric order parameters that respect the 
sympleetic nature of the phase space by considering the Lie algebra of generating functions of canonical 
transformations (symplectomorphisms) . 

The first step is to establish how a geometric quantity evolves in phase space, so that the sympleetic 
nature of its evolution is preserved. For this, we regard the action of the sympleetic algebra as an action 
of the generating functions h of k, rather than vector fields. The action is formally expressed as k ft, = 
£x^ K. The dual operation of the action which we denote by -k is then defined as ( k^Cj ^ ) = ( , ^ -^x^ C ) • 
Here Xh{q,p) is the Hamiltonian vector field generated by a Hamiltonian function h{q,p) through the 
definition _l w := dh. Notice that the star operation takes values in the space J-* of phase space 
densities kt^C G J-* . In the particular case of interest here, k is the phase space density k = f dqAdp and 
C = a function on phase space. In this case, the star operation is simply minus the canonical Poisson 
bracket, n-k g = [f^g\dq/\ dp. 

We shall first employ these considerations to find the purely dissipative part of the kinetic equation for 
a particle density on phase space. We choose variations of the form 5f = —£ Xh(4>) = ~ [/^(/): ^i'f)] 
with h{(j)) = {(j)-kf)^ = [(/),/] where ( • )" in (/Tk-^)' transforms a phase space density into a scalar function. 
The operation ( • will be understood in the pairing below. We then follow the steps: 
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Therefore, a functional F{f) satisfies the following evolution equation in bracket notation. 
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(3) 



The mobility /i and dissipation energy functional E appearing in (3) are modeling choices and must 
be selected based on the additional input from physics. The bracket (3) reduces to Kandrup's dissipative 
bracket for the modeling choice of = af with some a > [13]. The dissipation energy E in Kandrup's 
paper was taken to be the Vlasov Hamiltonian (see below), but in our approach it also may be taken 
as a modeling choice. This extra freedom allows for more fiexible interpretation and treatment of the 
dissipation process. 

Proposition 2.1 There exist choices of mobility ^[/] for which the bracket (3) dissipates energy E. 
Proof. The dissipative bracket in equation (3) yields E = {{E, E}} which is negative when /.j[./] is 
chosen appropriately. For example, /i[/] = /Af [/], where M[f] > is a non-negative scalar functional of 
/. (That is, M[f] is a number.) 

Remark 1 The dissipative bracket (3) satisfies the Leibnitz rule for the derivative of a product of func- 
tionals. In addition, it allows one to reformulate the equation (2) in terms of flow on a Riemannian 
manifold with a metric defined through the dissipation bracket, as discussed in more detail in [11]. 
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Proposition 2.2 (Casimir functionals) For an arbitrary smooth function (f> the functional C$ 
/$(/) dq A dp is preserved for any energy functional E. 
Proof. It suffices to cafcufate tlic bracket 
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(4) 



Corollary 2.3 The entropy functional S = // log / is preserved for any energy functional E. 
Remark 2 The existence of Casimirs and the corresponding preservation of any entropy defined solely 
in terms of f arises because the dissipative bracket (3) generates coadjoint motion, which is reversible. 
This property is shared with Kandrup's bracket, which is recovered for ^(/) = af for constant a > 0. 



3. Dissipation for kinetic moments: the Kupershmidt-Manin bracket and Darcy's law 



In this section we show how Eq. (2) leads very naturally to a nonlocal form of Darcy's law. In order 
to show how this equation is recovered, we introduce the Lie-Poisson structure for kinetic moments, also 
known as Kupershmidt-Manin structure [15]. We proceed by considering a one-dimensional configuration 
space; an extension to higher dimensions would also be possible by considering the treatment in [6]. 

As a general result [5,6,7], the equations for the moments of the Vlasov equation are a Lie-Poisson 
system. The n-th moment is defined as 



An{q) ■■= I p''f{q,p)dp. 
regulat 

{f,g}^Ia 



and the dynamics of these quantities is regulated by the Kupershmidt-Manin structure 

' 5F_ 5G_ 

5 An 5 Am 

where summation over repeated indices is omitted and the Lie bracket |-, ■] is defined as 

[am, /3„1 = n[3n{q)a'm{q)-mam{q)(3n{q) =■ ^^a^Pn 
The moment equations are 
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where the ad* operator is defined by (ad^^ Ak, otk-n+i ) {Ak, ad^^ otk-n+i )■ 

At this point one can consider the following Lie algebra action on Vlasov densities 

Pn f ■■= "^^-Yp.^„„ /=[./, P"Pn] (no sum) 

which is obviously given by the action of the Hamiltonian function h{q,p) = p"'Pn{q). Now, the dual 
action is given by 

{f^ng, Pn) ■■= (/, Png) = ( A 5 , (?) ) = (^J {f,g}p''dp, p, 
and the dissipative bracket for the moments (3) is written in this notation as 



{{E,F}} = -/jp" 



SE 



dp, p" 
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where we have substituted SE/5f = p'^Pk and SF/6f = p"''am and fis{q) Jp'^f^if] dp- 
Thus the purely dissipative moment equations are 

A„ = ad*^A,„+„_i with 7„i := (ad^^ Jlk+m-i)^ 



If we now write the equation for p :— Aq and consider only 70 and 71, we recover the following form of 
Darcy's law 



where we have chosen E = E[p] and JIq = p[p], so that 71 ~ po dqPo- 
3.1. Special cases 

Two interesting cases may be considered at this point. In the first case one makes Kandrup's choice in 
(1) for the mobility at the kinetic level p,[f] = f, so that Darcy's law is written as 



Kandrup's case applies to the dissipatively induced instability of galactic dynamics [13]. The previous 
equation is the Darcy law description of this type of instability. In the second case, one considers the 
mobility p[p] as a functional of p (a number), leading to the equation 



which leads to the classic energy dissipation equation, dE/dt = — (MPl^ffP)- 
3.2. Summary 

This section has provided a consistent derivation of Darcy's law from first principles in kinetic theory, 
obtained by inserting dissipative terms into the Vlasov equation which respect the geometric nature of 
the system. This form of the Darcy's law has been studied and analyzed in [9,10], where it has been shown 
to possess emergent singular solutions (clumpons), which form spontaneously and collapse together in a 
finite time, from any smooth confined initial condition. 



4. A dissipative Vlasov equation 

The discussion from the previous sections produces an interesting opportunity for the addition of dis- 
sipation to kinetic equations. This opportunity arises from noticing that the dissipative bracket derived 
here could just as well be used with any type of evolution operator. In particular, we may consider intro- 
ducing our bracket to modify Hamiltonian dynamics as in the approach by Kaufman [14] and Morrisson 
[16]. In particular, the dissipated energy may naturally be associated with the Hamiltonian arising from 
the corresponding Lie-Poisson theory for the evolution of a particle distribution function /. Therefore, 
we write the total dynamics generated by any functional F{f) as F — {F,H} + {{F, E}} where {• , ■} 
represents the Hamiltonian part of the dynamics. This gives the dissipative Vlasov equation of the form (2) 
with E = H, where H{f) is the Vlasov Hamiltonian. We illustrate these ideas by computing the singular 
(measure- valued) solution of equation (2), which represents the reversible motion of a single particle. 
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Theorem 4.1 Taking fi{f) to be an arbitrary function of the smoothed distribution f ~ K * f for some 

kernel K allows for single particle solutions f — 'Yld^i'^i^i.l ~ Qi{i))^{p ~ ^j(O)- The single particle 
dynamics is governed by canonical equations with Hamiltonian given by 

'SH 



{q,p) = {Q.{t),P.{t)) 



Proof. Let us write the equation of motion (2) in the following compact form 

with 



Sf 



and substitute the single particle solution ansatz f{q,p,t) ~ ^ Qi{t))^{p — Piit))- Now take 

the pairing with a phase space function (j) and write ( ^, /) = —([ 0, 7i ],/ ). Evaluating on the delta 
functions proves the theorem. 

Remark 3 The quantity —[ii{f),6H/Sf] plays the role of a Hamiltonian for the advective dissipation 
process by coadjoint motion. This Hamiltonian is constructed from the momentum map J defined by the 
* operation (Poisson bracket). That is, Jh{f,g) {g,-£x^f) = {g,[h,f]) = {h,[f,g]) = {h,f^g). 



5. Discussion and Conclusions 



This paper has developed a new symplcctic variational approach for modeling dissipation in kinetic 
equations based upon a double bracket structure in phase space. We focused our approach on the Vlasov 
example and found that the Vlasov case allows singic-particlc solutions, provided the mobility in the 
dissipation is a functional of the phase space distribution function. Moreover, we have shown how this 
approach recovers a nonlocal form of Darcy's law by using the Kupershmidt-Manin structure for kinetic 
moments. In general, it is also possible to extend our theory to the evolution of an arbitrary geometric 
quantity defined on any smooth manifold [12]. For example, the restriction of the geometric formalism for 
symplectic motion considered here to cotangent lifts of diffeomorphisms recovers the corresponding results 
for fluid momentum. One may also extend the present phase space treatment to include an additional set 
of dimensions corresponding to statistical internal degrees of freedom (order parameters, or orientation 
dependence) carried by the microscopic particles, rather than requiring them to be point particles. This 
is a standard approach in condensed matter theory, for example in liquid crystals, see, e.g., [2,3]. 

Remark 4 Being a special case of the dissipative bracket (3) presented in this paper, Kandrup's dou- 
ble bracket in (1) also possesses the Casimirs found in Proposition 2.2. However, the evolution under 
Kandrup 's double bracket does not allow single particle solutions. 

Remark 5 Had we chosen variations of the form Sf = —£ Xh(<l>) / = ^ [/; ^ ('/')] with h{(j)) = fJ-{f) * (/) 
— [/^(/) I 0] '^'^'^ followed the same steps as those in deriving (3), we would have obtained a different 
dissipative double bracket. It would have the same form as (3), but with <-> / switched in the 

corresponding entries. These two choices have different thermodynamic implications. In particular, the 
calculation in the proof of Proposition 2.3 would give entropy dynamics of the form 



5f 



M(/),log7 



/ - 

5f 



^0. 



For entropy increase, this alternative variational approach would require and E{f) to satisfy an ad- 

ditional condition (e.g., lJ-{f)/f and S E / S f functionally related). However, the Vlasov dissipation induced 
in this case would not allow the reversible single-particle solutions, because of the loss of information 
associated with entropy increase. 
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